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ABSTRACT 

The  elastic  displacement  field  of  moving  edge  dislocations  in  anisotropic 
body-centered-cubic  and  face-centered-cubic  crystals  is  found.  From  the 
elastic  displacements  the  shear  stress  on  the  dislocation  slip  plane  is  de¬ 
termined.  The  anomalous  velocity  range  in  which  edge  dislocations  of  like 
sign  attract  one  another  has  been  calculated  for  a  number  of  metals  and  ionic 
crystals.  It  is  found  that  anisotropy  does  not  appreciably  expand  the  anoma¬ 
lous  range. 

The  problem  of  dislocation  moving  on  the  interface  separating  media  of 
different  elastic  properties  has  been  considered.  The  anomalous  velocity 
range  may  or  may  not  exist,  depending  on  the  values  of  the  elastic  constants 
in  the  two  media.  The  dislocation  self-energy  is  infinite  at  the  slowest 
sound  velocity.  Supersonic  dislocation  behavior  is  qualitatively  described. 
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CHAPTER  I.  INTRODOCTION 


This  filial  report  describee  the  analysis  and  oaleulatlone 
on  fast  iBOvlng  dislocations  that  have  been  uislertaken  under 
this  contract.  Chief  interest  In  this  field  has  been  to  stu^ 
the  anoBalous  velocity  region  in  which  edge  dislocations  of 
like  sign  attract  rather  than  repel  one  another.  This  velocity 
region  occurs  from  the  Rayleigh  wave  velocity  (approximately 
0*9  times  the  slowest  sound  velocity)  to  the  velocity  at 
vMch  the  dislocation  self  energy  is  infinite.  In  the  anomalous 
velocity  region  dislocation  coalescence  can  take  place.  This 
coalescence  could  lead  to  crack  formation  and  fracture. 

Chapter  III  and  Chapter  IV  describe  the  results  found 
for  fast  moving  dislocations  moving  in  anisotropic  body-centered- 
cubic  crystals. 

Chapters  V  end  VI  present  similar  work  on  dislocations 
moving  in  face-oentered-oublc  crystals.  Finally,  in  Chapter  VII 
the  results  are  given  of  work  done  on  the  problem  of  dis¬ 
locations  isovlng  on  the  interface  separating  material  of 
different  elastic  properties. 

All  of  the  work  described  in  this  final  report  has  been 
written  for  Journal  articles,  and  has  appeared  or  will  appear 
in  the  llterat\ire  (VIeertman,  1962a,  1962b,  1962o}  Cotner  and 
Veertman  1962a,  1962b;  Van  Hull  and  Ueertman,  1962).  * 


*  References  are  listed  at  the  end  of  this  report. 
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Hore  oonplete  details  of  the  analTalB  and  caloulatiooe  are 
given  in  these  artioles  as  well  as  in  the  Status  Reports  1 
throng  6*  ** 


a 


**  Status  Reports  1  throui^  6  are  available  trcm  iFSMS 


CHAPTER  II  SUMMARY 

The  main  Interest  In  tMs  work  was  to  see  If  the  anomalous 
velocity  region  Is  larger  or  smaller  In  anisotropic  crystals  than 
It  Is  In  Isotropic  materials.  Erom  the  results  of  numerical  cal¬ 
culations  which  were  made  for  a  number  of  metals  and  ionic  crysteilB 
It  was  found  that  the  anomalous  region  Is  usimlly  smaller  In 
suolsotrpplc  crystals  than  In  isotropic  materials.  The  anomalous 
region  was  never  fouxti  to  he  appreciably  larger  than  that  of  Iso¬ 
tropic  crystals  althouj^  (Teutonlco,  1962a,  1962b,  1962c,)  It  is 
known  theoretically  that  It  Is  possible  for  the  regl(»i  to  be  very 
extended.  These  resxalts  are  similar  to  those  reeushed  by  Teutonlco 
(1962a,  1962b,  1962c)  through  numerlcsLl  calculations  concerning 
dislocations  different  from  those  considered  here. 

Since  dislocation  damping  forces  are  large  at  ordinary 
teiQ>eratures,  It  would  not  appear  likely  that  dislocations  would 
move  at  velocities  lying  In  the  anomalous  range.  It  may  be 
possible  to  have  dislocations  moving  In  the  anomalous  region  at 
very  low  temperatures,  where  dislocation  daniping  forces  become 
small,  or  at  very  high  stresses  such  as  are  encountered  In  shock 
loading  experiments.  If  dislocations  did  move  at  these  velocities, 
they  could  contribute  to  fracture  phenomena. 
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CHAPTER  in.  East  Movlog  Edge  Dislocations  on  tbo  (110) 

Plane  in  Anisotropic  Body-centered>ciibic 
Crystals. 

In  this  section  the  problem  of  an  edge  dislocation  moving 
uniformly  in  a  h.c.c.  lattice  is  analysed..  A  dislocation  is  con> 
sidered  which  lies  in  a  (110)  plane  parallel  to  a  < 112 > direction 
and  has  a  Burgers  vector  in  a  <111>  direction.  Edge  dislocations 
of  this  type  are  comnon  in  h.c.c.  crystals. 

It  has  been  pointed  out  (Weertnan  19^0,  1961)  that  edge 
dislocations  in  an  isotropic  medium  show  unusual  behavior 
idien  they  move  at  high  velocity.  Above  the  Raylei^  wave 
velocity  (approximately  0.9  times  the  transverse  soui^  velocity) 
the  stress  on  the  slip  plane  of  a  moving  edge  dislocation  changes 
sign,  and  edge  dislocations  of  like  sign  actually  attract 
rather  than  repel  one  another.  This  anomalous  behavior 
obviously  may  have  importance  in  fracture  phenomena  since  a 
coalescence  of  edge  dislocations  can  lead  to  crack  formation. 
Chief  interest  in  this  section  is  to  investigate  the  effect 
that  anisotropy  has  in  extending  or  contracting  the  velocity 
range  in  tdilch  the  anomalous  dislocation  behavior  occurs. 

The  problem  of  a  moving  edge  dislocation  was  first 
considered  by  Eshelby  (1949)*  ^dio  found  the  solution  for  the 
elastic  displacements  in  an  isotropic  medium.  Another  method 
of  obtaining  Eshelby 's  solution  was  given  later  by  Radok  (1956). 
BuUough  and  Bllbry  (1954)  next  considered  edge  dislocations 
moving  in  anisotropic  crystals.  They  found  the  solution  for 
the  elastic  displacements  for  the  case  in  which  the  problem 
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can  be  considered  as  one  of  plane  strain.  The  dislocation 
cannot  be  analyzed  as  a  problem  in  plane  strain»  and  a  more 
general  solution  for  the  elastic  displacements  must  be  found. 
Forttinately,  it  is  quite  obvious  from.  Bullough  and  Bllbgr's 
work  what  the  form  of  the  more  general  solution  must  be. 

THE»Kf 

Since  the  dislocation  under  consideration  lies  parallel 
or  peipendicular  to  the  crystal  directions  <111>»  <110>  and 
<112>,  it  is  oonnrenient  to  Introduce  a  coordinate  system  \dx>8e 
axes  run  along  these  directions.  Therefore,  we  shcLLl  adopt 
the  right-hand  coordinate  system  in  vdilch  the  positive  x  axis 
is  the  <110>  direction,  the  y  axis  is  the  <112>  ’•action, 
and  the  z  axis  is  the  <111>  direction. 


When  the  coordinate  axes  are  chosen  parallel  to  the 
three  <100>  directions,  a  cubic  crystal  has  three  independent 
elastic  constants:  c^,  ^  coordinate  system 

rotated  from  the  cube  axes  the  elastic  constants  are  given  by 
other  quantities.  Heannon(l956,  193?)  has  derived  e3q>licit 
formulas  for  obtaining  the  elastic  constants  in  any  rotated 
coordinate  system.  The  foUowizig  c  ustants  are  obtained  from 
Heaxmon's  formulas  for  the  coordinate  system  adopted  in  the 
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It  is  the  appearance  of  t  <v>n8tant  Cy^  In  the  nsatrix 
(l.l)  tMoh  makes  It  impossible  to  apply  BuUoujii  and  Billagr's 
solution  to  the  dislocation  tinder  study. 

The  equations  of  dynamic  equilibrium,  tdilch  must  bee 


satisfied  in  any  problem  in  elasticity,  are 


ijkl  2 


(sunned  on  J,k,  l,),..(l.3) 


®*k 


tdiere  p  is  the  density  of  the  crystcdj  t  is  time;  stanil.e  for 
x,y  or  s;  u^  stands  for  the  elastic  dlsplaoements  u,  r  and.  w  in  the 
x,y  and  s  directions  respectively;  and  are  the  ela»~fclo 

constants  in  the  unoontraoted  notation. 

Since  the  dislocation  being  considered  lies  parallel!  to 
the  y  axis,  the  derivatives  with  respect  to  y  can  be  set  ecjual 
to  sero;  otherwise  the  elastic  displacements  would  be  a  fnanotion 
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of  distance  along  the  dislocation  line.  Another  slsqplification 
in  the  equilibrium  equations  can  be  made  for  a  dislocation 
moving  at  a  uniform  velocity  V  in  the  z  direction.  The  solution 
of  any  of  the  three  elastic  displacements  u^  will  be  of  the 
foxm  u^  -  (x,s-Vt).  Hence  in  eqns.  (1.3)  the  time  derivative 

a^dt^  can  be  replaced  by  (d^/dz^),  Equations  (1.3)  are 
now  reduced  to  the  following  sett 
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(I.4a) 


(I.4b) 

(I.4o) 


In  contrast  to  the  case  of  an  edge  dislocation  moving  in  an 
isotropic  medium  or  to  the  problem  considered  by  BuUou^  and 
Billy)  it  is  Inqpossible  to  eliminate  one  of  these  equatioxis. 

All  three  elastio  displacements  are  involved  in  the  dlaplaoement 
solution. 

It  is  obvious  from  the  results  of  Bullou{^  and  Bllby 
that)  for  the  dislocation  being  considered)  the  solution  of 
elastic  displacements  which  will  satisfy  eqns.  (1.4)  very  probably 
has  the  fomt 

u  [A^log(8'^  +  r^x^)  +  k^g(z'^  +  Y2^x^)  +  A2log(a'^+Y3^x^)]  (I.5a) 
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(I.5b) 


V  » 


+  Bgtan”  ^  +  B^tan  ] 


2ff 


[C^tan“^IlJ  +C2taii“^I^  +  C^tan”^  ] . (I.5e) 


In  this  equation  s'  =  a  -  Vt,  b  is  the  length  of  the  Kirgers 
vector,  and  Bj^,  0^^,  Y^*  ®tc»  constants. 

To  evaluate  eqns.  (1.5)  idiloh  contain  12  unknown  constants, 
it  is  necessary  to  find  12  Independent  equations  in  these  constants. 
Nine  such  equations  are  obtained  )gr  setting  e^ns.  (1.5)  into 
eqzis.  (1.4).  For  exaople,  if  eqns.  (1.5)  are  placed  in  eqn. 

(I.4a),  the  following  equation  resultsi 
3  r 
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Since  this  equation  aust  bold  for  all  values  of  s'  and  x,  the 
bracket  term  containing  A^,  etc.,  must  equal  sero  for 
each  value  of  1.  Thus  from  eqn.  (1.6)  one  obtains  three 
eqiiatlons  in  the  unknown  constants.  Similar  substitutions 
into  eqns.  (I.4b)  and  (I.4o)  produce  six  more  equations.  The 
nine  equations  so  obtained  aret 


<"u'  -  -  '’’^>*1  *  n<“i3'  °  ^ 

Ti(cjj'  t  ■  ‘°33‘  “  ***’'’’' 
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(idiere  i  =  1,  2  or  3). 
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Since  the  right-hand  side  of  each  of  eqns*  (1*7)  is 
sero,  the  determinant  made  up  of  the  ooefflelents  of  etc, 

aunt  equal  aero  to  permit  all  the  equations  to  be  valid 
slmultaneouslyt  In  other  wordsi  must  be  one  of  the  three 
roots  of  the  oublo  equation 


•“2  ■"2 
(Y^  corresponds  to  the  minus  sign  and  to  the  plus  sign.) 

Equation  (1.8)  detamlnes  3  of  the  12  unknown  constants. 


EVALDATION  OF  THE  CONSTAHTS  A^^,  and  C^. 

The  constants  A^,  oan  be  evaluated  If  all  three  roots 

2 

Y-|^  are  positive  real  numbers  or  If  one  of  these  roots  is  a  posltlve.- 
real  number  and  the  other  two  are  oonqplex  conjugates.  Either  of 
these  two  possibilities  will  be  met  in  the  velocity  range  from 
V  -  0  to  the  velocity  at  vdilch  the  dislocation  self^onergy  becomes 
infinite. 
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2  2 

When  two  of  the  xx>ot8y  eay  Y2  f  oosqplex 

conjugates  (  and  hence  Yi  snd  are  likewise  cooqfilex  conjugates) 
it  is  only  neoessazjf.  that  the  pairs  of  constants  and  A^, 

«.d  Oj  Md  Oj,  ia«  b.  o<»pl«  eonjugate  nuiO^r.  to 
ensure  that  the  elastic  displacements  given  fay  eqn*  (1*5)  are 
real  numbers. 

To  evaluate  the  nine  other  tmknown  constants  appearing  in 
eqna.(l.5)B(l.6)  and  (1.7) >  it  is  necessary  to  obtain  three 
additional  independent  equations  containing  these  constants.  So 
far  the  condition  that  the  disj^cements  must  describe  a  dislocation 
has  not  been  used.  That  is,  if  a  circuit  is  made  around  the 
dislocation  line,  the  net  displacement  must  equal  the  Burgers 
vector.  In  eqns.  (1.5)  the  log  terms  return  to  their  original 
values  after  a  complete  elroxilt  is  made  around  the  dislocation 
line.  The  arc  tan  terms,  however,  change  by  a  factor  2it.  The 
Burgers  vector  of  the  dislocation  being  considered  lies  parallel 
to  the  z  direction.  Hence  it  can  be  seen  from  eqns.  (1.5)  that 
the  following  tvK>  equations  must  be  satisfied 

+  C2  +  ^  j.  .  .  .  .  (1.9) 

and 

Bj^  +  Bg  +  B^  =  0  .  .  .  .  (I.IO) 

Only  one  more  equation  is  needed.  It  is  obtained  from  the 
condition  that  point  forces  must  not  exist  at  the  core  of  the 
dislocation.  To  eliminate  the  point  forces  it  is  first  necessary 
to  determine  the  stresses  produced  by  the  elastic  displacements 
of  eqns.  (1.5).  The  following  set  of  equations  gives  the  stresses 
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in  the  coordinate  aysbem  being  used  if  all  derivatirea  with 
respect  to  y  are  aerot 
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To  ensure  that  no  point  force  acts  at  the  core  of  the  dislocation 
in  a  dlraetlon  perpendicular  to  the  slip  plana »  it  Is  only  necessary 
that  f  d  *0  •••  ••••  •  (1*12) 


/ 
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If  this  Integrel  were  not  equal  to  aero  a  point  force  would  have 
to  be  applied  to  the  core  of  the  dislocations.  The  magnitude  of 
this  force  would  be  eq\ial  to  the  integral.  When  eqns.  (1.3)  are 
substituted  into  eqn.  (I.U),  one  obtains  the  following  expression 
for  d. 
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Equation  (1.14)  is  the  result  of  substituting  this  e:q[>ression  into 
eqn.  (1.12)  (and  also  of  using  eqns.  (1.9)  and  (I.IO)  )t 
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Me  now  have  a  sufficient  number  of  equations  to  detexmlne  and 

C^.  Equations  (1.7) »  (I.9)i(l.l0)  and  (1.14)  can  be  rearranged  Into 
the  following  three  sets  of  equations^  from  which  it  Is  possible  to 
deteniilne  all  the  constants  A^»  and  The  sumnation  In  these 


equations  is  frcan  1  ^  x  to  1  =  3. 
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...(I.17a) 
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Z  ("44'  -  ZVi^ 
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The  solutions  of  eqns.  (I.15)b  (I.I6)  and(l.l7)  arei 
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(I.lda) 


®1  “  «(Y2^-Y3^)(b3Yi^  +  b2Yi%S^^» 


Bg  =  *(Y3^  -  Yi^)(b2Y2^  +  b2Yi\2%^)»  •• 

Bj  =  .(t/  -  r/XbjT^^*  *  b^nSV),  .. 

\diere 

»2  =IVl"- 

‘'3  =IW' 

“  Y3^)(83Yi^  +  ~  (l«19a) 

°2  “  “^Y3^  -  ■*■  “  '^Z'^\  "  Y2^Y3^)»  (I.i9b) 

S  *  ■*■  *2^1^^2%^  "  "  ^3^^2^^»  (I«19c) 

\didre 

«2 

«3 

Aj^  =  *Yi'’^(Y2^  -  Y3^)(a3Yi^  +  Y3^))»  (l»20a) 

^2  =  •Y2“^^'^3^  “  ^  ■*■  ■^3^^^* 

Aj  =  «Y3"^(Yi^  “  ^2^^ ^*3^3^  '*'  ^  (l.20o) 

where 

^  “  °13^  /  ®n' 

*^2 

*3  “liVl^* 
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A  eonaideratlon  of  these  eqiiatlons  will  show  that  it  is  indeed  true 

that  vdien  Y2  ooiqilex  oonjugates,  the  constants  and  A^y 

etc.,  also  are  ooiq>lex  conjugates.  Equations  (1. 18),  (1.19)  and  (1.20), 

along  with  eqn.  (I. 8),  give  the  general  solution  of  our  problem  in  the 

velocity  range  from  V  =  0  to  the  velocity  at  which  the  dislocation 

self-energy  becomes  infinite. 

SHEAR  STRESS  ON  THE  DISLXATION  SLIP  PLANE 

Of  all  the  stresses  around  an  edge  dislocation,  the  stress 

\dtloh  is  of  primary  interest  is  that  stress  which  produces  a  force 

on  another  edge  dislocation  with  a  parallel  Burgers  vector.  This 

force  is  <s  b,  tdiere  b  is  the  Burgers  vector.  Hence  0  is  the  shear 
xa  xz 

stress  of  greatest  interest  to  us.  This  stress  can  be  found  Iqr 
substituting  eqns.  (1.5)  into  eqn*  (l.U).  One  obtaims 


xs 


“  2n  ^1  ^2  2  2  ^  (I.21a) 

*  +  Y^  * 


The  shear  stress  on  the  slip  plane  of  the  dislocations  Itself  (3PO)  is 

^'’xa^xO  "  2^^.  L  ^®44  *i  *  ^i  °14  ®1  ^44  °i  ^  (I.21b) 


In  the  case  of  a  stationary  or  slowly  moving  dislocation  this  stress 
can  be  expected  to  be  positive  for  a  positive  edge  dislocation. 

Slowly  moving  dislocations  of  like  sign  and  on  the  same  slip  plane 
repel  one  another.  It  is  known  (Usertman  1961)  that  an  edge 
dislocation  moving  in  an  isotropic  medium  e^^rlenoes  a  shear  stress 
on  its  slip  plane  idilch  decreases  with  increasing  dislocation  velocity, 
and  which  actually  can  change  sign.  In  the  velocity  region  where  the 
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sign  of  the  shear  stress  has  been  reversed,  dislocations  of  nVo 
sign  attract  rather  than  repel  one  another.  The  velocity  separating 
this  region  of  abnomal  behavior  from  the  normal  region  is  that  velocity 
at  vihioh  the  shear  stress  is  exactly  equal  to  zero.  This  velocity 
turns  out  to  be  the  Raylelf^  vave  velocity. 

To  ensure  the  existence  of  an  abnormal  region  in  \dileh 

dislocations  of  like  sign  attract  one  another  it  is  necessary  only 

that  eqn.(l.21b)  goes  to  zero  at  a  velocity  lower  than  the  velocity 

2 

at  vdxloh  one  of  the  roots  Is  equal  to  zero  or  at  \diieh  two  of 
the  roots  beocmis  negative.  The  velocity  at  vdilch  eqn.  (1. 21b)  is 
equal  to  zero  is  the  velocity  for  \Mch 

^1®14^®1  ^i®44^*'i  /  “  0  •  •  •  •  (I»22a) 

The  velocity  satisfying  this  equation  is  the  velocity  at  idilch  a 
generalized  Rayljij^  wave  >K>uld  travel  on  a  (llO)  plane  in  a  <111> 
dlrestion.  Since  tho  shear  stress  o  ,  in  general,  is  not  equal 
to  zero  in  the  dislocation  slip  plane,  an  actual  surface  wave  cannot 
propagate  on  a  (110)  plane  in  a  <111>  direction  in  an  anisotropic 
crystal. 

Equation  (I. 21a)  can  be  slnqilified  somewhat  throu^  the  use 
of  eqn.  (I. 7b).  If  the  latter  equation  is  divided  through  by  snd 
then  Bunmed  from  1  =  1  to  i  =  3  and  the  result  substituted  Into 
eqn.  (I. 21a),  one  obtains 

?  [“13'*!  "  ^^33  "  Vl'^]  =  °  . 

This  equation  has  the  advantage  that  all  the  constants  have  been 
eliminated.  If  eqns.  (1.19)  and(l.20)  are  substituted  into  eqn.  (I. 22b) 
the  following  equation  results! 
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where 


I'l'l'j  *  TjY3  *  T3T1) 

*  *2'f-T2T3<ri  +  Yj  +  T3)  -  13  “  0  •  ••  •••  •••  (l.22o) 

0 

-  ^*33'  “  -  ^^33  -  P^>» 


“U 


^2  “  °13'‘*2  ■  *°33*  ■  f''^)*2  "  < 


'U  -  o’" 


°13 


'3  “  «13'‘3  -  <'3l'  -  0*^)*3  =  -  ^  <°U'  *  V  -0  O’*). 


"u 


+  2(0^2'  -  pV^)  -12^+ 
®11 


(033'  -  pV^)  +  (c,,'  /  o^/)(o,/  -  O,,') 


13  '  ''U  13 


®U  66  ■  ®U 


"^/^®66*°13^  ®66*®44^  “  ^®U  "  ®66^^®33'  “ 


^11 


The  principal  goal  of  this  section  was  the  drivation  of  eqn* 

(I.22o).  MLth  this  equation,  one  can  inrestigate  the  effect  of 

anisotropy  on  the  extent  of  the  velocity  region  in  which  anmalous 

dislocation  behaviour  occurs*  For  an  anomalous  velocity  region 

to  exist  it  is  necessary,  of  course,  that  the  velocity  \Moh 

satisfies  eqn.  (1.22)  is  smaller  than  the  velocity  at  tMch  the 
2 

roots  first  become  aero  or  negative.  The  dislocation  energy 
is  inflxilte  at  this  latter  velocity  and  hence  it  is  the  limiting 
velocity  of  dislocation  motion  in  noznal  circumstances.  In  Chapter 
IVresults  of  numerical  calculations  using  Eqn.  (1.22)  are 
presented. 
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CHAPTER  IV.  Calculations  for  Body-CenteredTpublc  Crystals 

This  section  presents  the  results  of  a  numerical  calculation 
of  the  shear  stress  on  the  slip  plane  of  an  edge  dislocation 
moving  uniformly  in  a  body-centered  cubic  lattice.  The  dislocation 
considered  is  moving  on  a  (UO)  plane  in  a  <111>  direction. 

The  primary  concern  is  to  investigate  the  extent  of  the  anomalous 
velocity  region  in  idiloh  the  shear  stress  reverses  sign  and 
dislocations  of  like  sign  attract  one  another.  This  anomalous 
behavior  of  the  dislocations,  \dilch  leads  to  a  coalescence  of 
fast  moving  dislocations,  may  be  expected  to  be  of  ljqx>rtance  in 
fracture  phenomena.  Iron,  therefore,  is  logically  the  metal  to 
study. 

In  the  previous  section,  formulas  were  developed  iMch 

permit  the  calculation  of  the  shear  stress  on  the  slip  plane 

of  an  edge  dislocation  moving  on  a  (UO)  plane.  It  was  found 

that  the  shear  stress  d  acting  on  the  slip  plane  at  a  distance 

z'  from  tlie  center  of  a  dislocation  moving  with  the  velocity  V 

is  0  =  -  "^2  '*'3^  ^1^2  V3  '^3^1  ^ 

2Tt8'  (Il.l) 

+  f2YiY2Y3(Yi  +  Y2  +  Y3)  -  ^3  ■*-(Yi^Y2)  (Y2-^3)  (Y3ni) 

where  b  is  the  Burgers  vector  and 
/2 

fj  =  -  (033'  -  pV^) . (II*2a) 

®ll' 

^2  “  . (n.2b) 
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f3  =  —  ( V'  ^®13'  +  —  <“33'  - 


®11 


“n 


(C„'  -  pv")  .  !]tL  (C 


“33  -  ^  ■  — -r  “  *13' 

_ !ll _ 

®44  66  "  ®14  ^ 


“13 


/  ^“66^3^  ®66*®44*  “  "  ®66^^®33*“P^) 


“11 


2  2  2 

and  »  Y2  Y3  are  t)ia  three  roots  of  the  oublo  equation 
given  hgr  Eqn.  (1. 8). 

RESULTS  FOR  IRON 


In  Table  I.l  are  listed  values  of  these  elastic  constants 
for  iron  calculated  from  data  of  Seits  and  Read  given  in  Heanaon's 
review  article  (1946). 

These  values  were  used  to  calculate  these  roots  of  equation 
(1.3) «  which  are  tabulated  in  Table  11.2. 

The  limiting  dislocation  velocity  is  found  to  be  0.941o» 
\diere  c  =  *^44  Vp)  =  the  velocity  of  shear  waves  in  a  <111> 
direction  in  a  cubic  crystal. 

The  limiting  velocity  of  dislocation  motion  in  iron  is 
smaller  than  the  shoar  wave  velocity  in  the  direction  of 
dislocation  notion.  The  possibility  that  the  dislocation  motion 
may  be  limited  to  a  velocity  smaller  than  the  shear  wave  velocity 
first  was  pointed  out  by  L,  J,  Teutonioo  (I962ay  1962by  1962o). 

The  shear  stress  on  the  slip  plane,  plotted  in  Fig.II.l 
was  calculated  from  equation  (I.l)  by  using  the  values  of  the 


(I1.2c) 
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TABLE  I.l 


11  2 

Elastic  constants  of  iron  (in  units  of  10  dyVin  .) 


®11 

®12 

®U 

23.7 

14.1 

11.6 

°ll' 

®33* 

\2 

®13* 

°44* 

«66' 

30,5 

32.8 

11.8 

9.56 

7.06 

9.33 

3.40 
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0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.3 

0.82 

0.84 

0.86 

0.88 


4.77 

4.61 
4.46 
4.33 
4.16 
4.04 
3.90 
3.76 

3.62 
3.59 
3.57 
3.54 
3.53 


2.2 

Yi  BxAXz 

0.359  ±  0.2781 
0.325  t  0.2681 
0.277  t  0.2501 
0.226  ±  0.2311 
0.180  t  0.2131 
0.133  1  0.1911 
0.086  ±  0.1691 
0.038  ±  0.1351 
■0.008  ±  0.09251 
-0.020  ±  0.07971 
.0.029  ±  0.07101 
.0.038  t  0.05201 
-0.048  ±  0.032291 


TABLE  11.3 


0 

0.269 

0.1 

0.259 

0.2 

0.251 

0.3 

0.2a 

0.4 

0.229 

0.5 

0.214 

0.6 

0.190 

0.7 

0.160 

0.8 

0.0814 

0.82 

0.0568 

0.84 

0.00925 

0.86 

-0.0795 

0.88 

-0.315 
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PIG,  II.l.  Plot  of  shear  strasa  a  on  slip  plane  versus  pV  /o^* 


22 


TABLE  II. 4 

Calculated  Charaoterlatlo  Dislocation  Velocities*  * 


Material 

°i 

Ac 

Body-oentered-cubic  crystals 

U 

O.8O9 

0.803 

0.006 

1.042 

P  Brass 

0.824 

0.815 

0.009 

1.010 

Na 

0.869 

0.855 

0.014 

0.965 

K 

0.853 

0.842 

0.011 

0.905 

oFe 

0.941 

0.918 

0.023 

0.485 

W 

1.000 

0.913 

0.087 

0 

Mo 

1.000 

0.903 

0.097 

-0.115 

Csl 

1.000 

0.835 

0.165 

-Okli56 

CsBr 

1.000 

0.895 

0.105 

-0.167 

RbBr 

0.754 

0.710 

o.ou 

-0.439 

Rbl 

0.726 

0.693 

0.033 

*0*464 

oFe^®^ 

0.952 

0.927 

0.026 

0.439 

oFe^^^ 

0.949 

0.923 

0.024 

0.433 

Li^®^ 

0.809 

0.804 

0.005 

1.042 

♦Velocities  are  expressed  in  units  of  (o^Vp)'*’^  =  transverse 

sound  velocity  in  the<lll>  direction.  Elastic  data  used  in 

calculations  are  those  listed  in  Huntington's  review  article  (1958) 

for  Li,  K,  Na  (data  of  Bender),  U,  Ho,  and  ^  brasss  (data  of  Artman 

and  Thoiiq>son);  in  Reluts^s  paper  (1961)  for  Csl,  CsBr,  RbBr  ,  and 

Rbl;  and  in  Heamon's  review  article  (1946)  for  oFe  (data  of  Belts  and  Read). 

(a)  Calculated  by  L,J.Teutonlco  (private  communication)  using  elastic 
data  listed  in  Heannon  (1946)  (data  of  Kimura  &.  Ohno) 

(b)  Calculated  by  L. J .Teutonlco  (private  communication) using 

elastic  data  of  Bayne  and  Chandrasekhar  (1961)  (footnote  cont.  on  next  page) 
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rootB  of  equation  (1.8).  The  actual  calculated  Talues  of  the 
shear  stress  are  listed  in  Table  II. 3.  It  can  be  seen  that  the 
stress  is  positive  at  low  velocities)  goes  throuf^  zero  and  then 
becomes  negative  at  &i^  velocities.  The  velocity  at  vMch  the 
stress  is  zero  is  the  Rayleigh  mve  velocity.  Its  value  is  O.9I80. 

The  result  that  the  anomalous  velocity  region  extends  from 
only  O.9I80  to  O.94I0  is  in  qualitative  agreement  with  calculations 
(tfeertman)  1962a)  made  for  dislocations  moving  in  slightly 
anisotropic  lattices.  From  these  calculations  It  was  found  that 
a  slight  anisotropy  (with  c^'  taken  to  be  a  positive  number) 
decreases  the  extent  of  the  anomalous  velocity  region.  An 
isotropic  crystal  with  the  sane  <^44 ratio  as  iron  has  a 
Rayleli^  wave  velocity  of  the  order  of  0.94c.  The  limiting  velocity 
of  dislocation  motion  ior  an  isotropic  crystal  is  always  o. 

Hence  the  extent  of  the  anomalous  velocity  range  for  edge  dislocations 
moving  on  the  (UO)  plane  in  iron  is  3  times  smaller  than  that 
of  an  Isotropic  crystal  with  cosy^arable  values  of  and 
Thus  it  is  more  difficult  in  iron  to  bring  dislocations  moving 
on  a  (no)  plane  into  a  velocity  range  \diere  dislocation  coalescence 
can  take  place. 

RESULTS  ON  OTHER  BODY-CENTERED-^BIC  CRYSTALS 

Table  II. 4  indicates  the  principal  results  on  calculation  on 
other  b.o.o.  crystals.  This  table  lists  the  velocity  c^^  at  which 
tne  self -energy  of  a  moving  dislocation  is  infinite)  the  velocity  o^, 


(footnote  continued  from  previous  page) 

(0)  Calculated  by  L.J.  Teutonico  with  elastic  data  listed  in 
Huntington  (1958). 
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(the  generalized  Rayleigh  wave  velocity)  at  idilch  the  shear  stress 
on  the  dislocation  slip  plane  Is  zero,  the  difference  between 
these  velocities  Ac  =  c.  -  c,  and  the  ratio  of  the  elastic  constant 
where  c^^^  and  c^*  are  elastic  constants  In  the  rotated 
coordinate  system.  The  ratio  is  a  measure  of  the  degree 

of  anisotropy  of  the  elastic  constants.  The  quantity  c^^*  is  equal 
to  zero  for  an  Isotropic  material.  Listed  In  Table  11.4  are  the 
results  of  calculations  by  Teutonico  on  alpha  iron  and  lithium. 

In  the  velocity  range  from  c^  to  c^  dislocations  on  the 
same  slip  plane  of  like  sign  attract  rather  than  repel  one  another. 

In  this  velocity  range  dislocation  behavior  Is  anomalous. 

From  calculations  of  the  effect  of  a  slight  anisotropy 
(Weertroan,  1962a)  on  the  extent  of  the  anomalous  velocity  range 
it  was  predicted  that  when  is  positive,  an  increase  in  the 
anisotropy  decreases  the  anomalous  velocity  range ,vrtierea8  if  c.j 
Is  negative,  the  anomalous  range  Increases.  An  inspection  of 
Table  11.4  reveals  that  crystals  with  a  positive  c^'  do  have  a 
smaller  Ac  than  ttmgsten,  which  is  an  almost  isotropic  material. 

The  crystals  Csl  and  CsBr,vdilch  have  a  negative  ,have  a  larger 
Ac  than  tungsten.  However,  RbBr  and  Rbl,  vdilch  are  more  anisotropic 
than  Csl  and  CsBr  and  >diich  also  have  a  negative  c^' ,  have  a  smaller 
Ac  than  tungsten.  It  is  clear  that  the  resrilts  of  the  sli^tly 
anisotropic  calculations  cannot  be  extrapolated  to  large  values 
of  anisotropy  and  that  there  is  a  limit  to  the  extent  to  vAilch 
anisotropy  can  widen  the  anomalous  velocity  range. 
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CHAPTER  V.  F&st  Moving  Edge  Dislocation  on  the  (ill)  Plana  in 
Anisotropic  Face-Centered>Cubic  Crystals 

This  section  considers  the  problem  of  edge  dislocations 
moving  uxilfoznly  in  anisotropic  faoe-centered-cubic  crystals. 

Considered  is  the  ordinary  edge  dislocation  of  face-centered-cublc 
crystals  lAileh  lie  in  a  (ill)  plane  parallel  to  a  <112>  direction 
and  have  a  Burgers  vector  in  a  <110>  direction.  In  the  following 
section  the  elastic  displacement  field  of  this  type  of  dislocation 
will  be  determined  and  also  the  shear  stress  on  the  dislocation 
slip  plane.  The  anomalous  velocity  range  can  then  be  detemlned 
from  this  shear  stress. 

ELASTIC  DISPLACfMEMT  FIELD 

The  coordinate  system  is  adopted  in  vdiloh  the  x  axis  is  paxmllel 
to  the  <110>  direction,  the  y  axis  is  parallel  to  the  <112>  direction, 
and  the  %  axis  is  parallel  to  the  <111>  direction.  The  elastic 
constants  in  this  coordinate  system  are  the  same  as  those  given  in 
Chapter  HI. 

Since  the  dislocation  being  considered  lies  parallel  to  the 
y  axis,  the  elastic  displacements  about  the  dislocation  must  be 
independent  of  y.  Since  a  dislocation  moving  uniformly  with  a 
velocity  V  in  the  x  direction  is  also  being  considered,  the  elastic 
displacements  will  be  a  function  of  x  -  Vt,  where  t  is  the  time. 

If  u,  V,  and  w  represent  the  elastic  displacements  in  the  x,  y, 
and  z  directions  respectively,  the  equations  of  dynamic  equilibrium 
reduce  to  the  following: 
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It  Is  obvious  from  the  results  of  BuUough  and  Bilb7  (1954) 
that  the  solution  of  the  elastic  displacements  vdiich  satisfies  Bqns* 
(lll.l)  is  very  likely  to  have  the  form: 


u 


V 


^  I  A^tan 

2ii  \ 

__  Bj^logCx'^  +  Y^^z^)  +  B2log(x'^  +  r2^*^)  +  B^log(x'^  + 

4ff 


.lYi* 


Agtan' 


-iV 


(III.2a) 


w  = 


4n 


jC^log(x'^  +  +  C2log(x'^  +  ^2^*^)  Cjlog(x'^  + 


where  x'  =  x  -  Vt,  b  is  the  length  of  the  Burgers  vector,  and 
B^,  Yj^i  etc*,  are  constants.  These  equations  contain  12  unknown 
constants  which  must  be  determined.  Equations  (III. 2)  are  similar 
to  a  set  which  give  the  elastic  displacements  about  an  edge 
dislocation  moving  on  a  (llO)  plane  in  a  body-oentered-«ublo  crystal. 
The  12  unknown  constants  can  be  evaluated  in  a  manner  similar  to  that 


(I11.2b) 

(lII.2o) 
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carried  out  in  the  body*^enterad->ouble  case.  This  evaluation  has 
been  carried  out  by  Usertman  (1962b). 

SHEAR  STRESS  ON  THE  DISLOCATION  SLIP  PLANE 

The  stress  vdilch  produces  a  force  on  another  parallel  edge 
dlslooatlon  with  a  parallel  Burgers  vector  Is  the  shear  stress 


8u  av  aw 


(III.3) 


vMch  with  appropriate  substitutions,  becomes 
b  ^  x‘ 

‘’xz  “  —  L  —5 - 5-5-^®U  ^1*1  ®U  ®i  *^1^  (III»4) 

2  1  X  +  Yi  a 

On  the  slip  plane  of  the  dislocation  Itself  (z  =  0)  this  equation 
reduces  to 


’xz  “ - L  ^®44*  Vl  ®U  ®i  *^1^ 

2irx'  1 


(ni.5a) 


which  also  can  be  written  as 


b  r 

®xa  “ -  E  ril^  "  \  ~  ®14*®1  ®13'°i  ^ 

2irx'  1 

If  the  values  of  A^,  and  (Weertman  1962b)  are  substituted  In 
}!iq.(lll.5b)  one  finds! 


(111.5b) 


2nx'  (Yi-^2^^^2'*^3^^^3'^1^ 


where 
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At  the  Rayleigh  viave  velocity  the  shear  stress  eqtials  aero. 

At  velocities  below  the  Raylel^i  wave  velocity  will  be  positive 

aiid  above  it  this  stress  will  be  negative.  In  the  velocity  range 

tAiere  a  is  negative,  dislocations  of  like  sign  attract  rather  than 
xz 

repel  one  another.  Thus,  the  velocity  at  \dilch  is  zero  is  of  great 

Interest;  the  behavior  of  dislocations  moving  slower  than  this 

velocity  la  normal  vdiereas  dislocations  moving  faster  than  this  velocity 

e^diibit  abnonool  behavior. 

LIMITING  VELCX3ITY  OF  DISLOCATION  MOTION 

As  the  velocity  of  a  dislocation  increases,  its  self-energy 

also  increases.  At  some  velocity  the  self-energy  will  become 

infinite,  and  this  velocity,  therefore,  sets  an  upper  limit  to 

dislocation  motion.  It  was  previously  noted  (Ueertman  1962a)  that 

the  limiting  dislocation  velocity  is  the  smallest  velocity  at  vdildh 

one  of  the  roots,  Yj  >  ^'Irst  becomes  zero  (if  all  roots  are  real 

nvimbers)  or  becomes  a  real  but  negative  nvimber  (if  two  of  the  roots 

are  complex  conjugates).  The  velocities  at  which  the  roots  become 

2 

zero  can  be  fovind  by  setting  Yj|^  equatl  to  zeiD  .  Thus  one  can 
obtain  the  following  equation: 

(06^'  -pv2)(c^'  -pv2)(ci3^'  -pV^)  = 

The  smallest  velocity  satisfying  this  equation  is  the  velocity  at 
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2 

^oh  one  of  the  roots  first  becomes  zero.  The  sazUest  veloolty 
satisfying  Eqn.(III.6)  Is  glTsn  by  the  quadratlo  equation 


{py2)2  .  pt2(o^'  +  o^^')  +  .  oj^'2  =  0  (ni.7.) 


vdiose  solution  Is 

\  *  =66'>  - 

Jz 


l=u'l 


(111.7b) 


If  Is  positive 


>  ^  this  equation  can  be  written  as 

'u'  / '46'  -  /V 


If  Is  negative  >  °66^^ 


The  velocities  s^  and  S2» 
direction  are  (Viateman  1959) 


'a  '  |°u  ■  '12'/^ 


+  /2C14'  (in.7d) 

of  the  two  shear  waves  In  <110> 


(Ill.Sa) 

(111.8b) 


In  terms  of  c 
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®66'»  ®u' 


these  eqviatlons  become 


“1  =  H'u  =  t‘'66' 

•2  =  |°u'  - =  ['66'  -y^  v'*/' r  ^  <“*•“> 


If  c^'  Is  positive,  S2  Is  the  slower  shear  wave  velocity,  and  If 
c^'  Is  negative,  s^^  Is  the  slower.  From  Eqns.  (III.7c)  and 
(III.7d),  It  can  be  seen  that  a  root  first  becomes  zero  at  the 
slower  of  the  two  shear  wave  velocities  in  the  <11Q>  direction.  I& 
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2 

the  slightly  anisotropic  case  the  roots  always  positive  mmbers 

idian  the  dislocation  velocity  is  less  than  either  s^  or  S2  hut  is 
not  too  close  to  zero.  Hence,  for  this  case  the  slower  of  the  two 
shear  wave  velocities  is  the  limiting  velocity  of  dislocation  siotion. 

In  the  general  anisotropic  case,  two  of  the  roots  can  be  complex 
conjugates  and  the  dislocation  limiting  velocity  will  he  the  smallest 
velocity  at  vdilch  either  coaQ>lex  conjugate  roots  turn  into  a 
negative  number,  or  a  root  becomes  equal  to  zero. 
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CHAPTER  VI.  Calculations  for  Faoe-^entered-^Jublo  Crystals 

In  Chapter  V  an  equation  was  derived  \diioh  gave  the  shear 
stress  acting  on  the  slip  plane  of  a  moving  edge  dislocation  in  an 
anisotropic  face-centered«oabic  lattice.  This  section  presents  the 
results  of  numerical  calculations  of  the  shear  stress  on  the 
slip  plans  using  that  equation.  The  metals  tdilch  have  been 
investigated  are  aluminum,  copper,  gold,  silver,  lead,  and  nickel. 

These  calculations  are  conqilementary  to  those  on  body-oentered-cubic 
crystals. (Chapter  IV) 

Table  IV.l,  lists  the  elastic  constant  data  used  in  the 
calculations.  These  data  are  taken  from  Huntington's  review 
article. (1958). 

In  Table  IV. 2  are  listed  calculated  values  of  the  shear 

stress  on  the  slip  plane  at  various  velocities.  These  shear  stresses 

are  plotted  as  a  function  of  dislocation  velocity  in  Figs.  IV.l  and 

IV.2.  The  most  striking  result  contained  in  these  figures  and  the 

table  is  the  extreme  narrowness  of  the  velocity  range  in  vnilch  the 

shear  stress  is  negative.  (It  is  in  this  range  that  dislocations 

of  like  sign  attract  one  another,  and  thus  dislocation  coalescence 

can  take  place).  Aluminum,  however,  is  an  exception;  here  there  is 

an  appreciable  velocity  region  in  idilch  the  shear  stress  is  negative. 

Aluminum  has  such  an  extended  anomalous  region  simply  because  it 

is  an  almost  isotropic  material.*  It  is  known  (Weertman,  1981)  that 

*Huntington  lists  cmother  set  of  elastic  constants  for  alumimim  idiich 
are  slightly  more  anisotropic  than  those  used  in  the  present  paper. 

We  carried  out  calculations  using  this  more  anisotropic  data  and 
concluded  that  the  anomalous  region  is  not  appreclaU.y  reduced  in  sin.  W2 
It  was  foxmd  that  c  =  0.944  and  v^  =  0.919(in  velocity  units  of  (V  /p)^  . 
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TABLE  IV.l 

U  2 

Elastic  constants  of  metals^  (In  units  of  10  dynes/cn  ) 


A1 

Cu 

Au 

Pb 

Ni 

Ag 

10.8 

16.8 

18.6 

4.66 

24.7 

12.4 

®12 

6.13 

12.1 

15.7 

3.97 

U.7 

9.34 

"44 

2.85 

7.54 

4.20 

1.44 

12.5 

4.61 

11.3 

22.0 

21.4 

5.73 

32.2 

15.5 

®33^ 

U.5 

23.8 

22.3 

6.09 

34.7 

16.5 

"U 

2.51 

4.08 

2.37 

0.73 

7.46 

2.56 

2.68 

5.81 

3.28 

1.08 

9.97 

3.58 

°12^ 

5.96 

10.4 

14.8 

3.56 

12.2 

8.31 

"13' 

5.79 

8.68 

13.9 

3.21 

9.72 

7,29 

"u 

0.24 

2.45 

1.30 

0.50 

3.54 

1.45 

a  Data  taken  from  Huntington  (1958) 
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table:  IV.  2 


Calculated  values  of  the  shecur  stress  on  the  slip  plane  at 
various  velocities. 

Nickel 


pV^ 

dx' 

0 

0.263 

0.1 

0.255 

0.2 

0.246 

0.3 

0.236 

0.4 

0.223 

0.5 

0.204 

0.6 

0.161 

0.64 

0.0883 

0.65 

0.0368 

0.66 

-0.123 

0.665 

M  CO 
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TABLE  17.2  (Continued) 


Lead 


pV^ 

ax' 

'-u' 

0 

0.324 

0.1 

0.317 

0.2 

0.310 

0.3 

0.301 

0.4 

0.286 

0.5 

0.186 

0.505 

0.119 

0.507 

o.oas 

0.5075 

-0.00340 

0.508 

-0.0721 

0.50928 

—  09 
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TABLE  IV*2  (Continued) 
Gold 


^44 

0 

0.330 

0.1 

0.325 

0.2 

0.315 

0.3 

0.305 

0.4 

0.301 

0.5 

0.271 

0.6 

0.U7 

0.608 

0.0225 

0.61275 

-  CO 

36 


TABLE  IV. 2  (Continued) 


Aluminum 


pV^ 

<jx' 

^44' 

0 

0.255 

0.1 

0.247 

0.2 

0,237 

0.3 

0.226 

0.4 

0.213 

0.5 

0.197 

0.6 

0.176 

0.7 

0.144 

0.8 

0.099 

0.87 

-0.0103 

0.9 

-0.156 

0.93303 

»  ca 
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TABI£  IV.2  (Continued) 


Copper 


pV^ 


V 

0 

0.313 

0.1 

0.312 

0.2 

0.306 

0.3 

0.300 

0.4 

0.291 

0.5 

0.261 

0.56 

0.161 

0.57 

0.093 

0.572 

-0.06370 

0.574 

-0.119 

0.57591 

.  «a 
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TABLE  IV.2  (Continued) 


Silver 


4 

ax! 

°44 

0 

0.316 

0.1 

0.310 

0.2 

0.304 

0.3 

0.296 

0.4 

0.286 

0.5 

0.268 

0.58 

0.212 

0.59 

0.133 

0.595 

0.0134 

0.59839 

.  ca 
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on  extensive  anomalous  velocity  region  exists  In  Isotropic  materials* 

A  slight  anisotropy  usually  decreases  (Weertman,  196^1^ 1962b) 
the  anomalous  velocity  region.  Therefore,  It  Is  not  surprising 
that  the  anomalous  region  Is  decreased  In  extent  by  a  large  anisotropy. 
A  reduction  was  found  In  the  case  of  alpha  Iron  (Gotner  &  Ueertman, 
1962a,  1962b)  \dilch  Is  a  strongly  anisotropic  body-centered-cublc 
metal.  However,  the  degree  of  decrease  in  the  strongly  anisotropic 
fee  metals  Is  much  larger  than  it  Is  In  alpha  Iron.  To  Illustrate 
this  point.  Table  IT. 3  lists  the  velocity  at  which  the  anomalous 
region  starts  (v^  =  the  Rayleigh  wave  velocity)  and  ends  (c  =  the 
velocity  at  vdilch  the  dislocation  energy  Is  Infinite). 

From  Table  IV. 3  It  Is  clear  that  to  bring  an  edge  dislocation 
Into  the  velocity  region  in  vdiich  dislocation  coalescence  occurs 
would  be  more  difficult  in  the  case  of  strongly  anisotropic  foe 
metals  than  of  strongly  anisotropic  bcc  metals.  It  would  be 
tempting  to  ascribe  the  greater  ease  of  fracture  of  bcc  metals 
to  this  difference*  However,  the  fact  that  aluminum  has  a  large 
anomalous  region  and  yet  appears  to  be  ductile  down  to  low 
temperatures,  presents  a  serious  difficulty  to  such  a  theory. 
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TABLE  IV. 3® 


Metal 

0 

aFe 

0.94 

0.92 

A1 

0.966 

0.933 

Cu 

0.759 

0.756 

Au 

0.783 

0.780 

Pb 

0.7U 

0.712 

Hi 

0.815 

0.806 

Ag 

0.774 

0.771 

1/2 


a  Velooitles  are  esqpressed  in  unite  of  (o^Vp) 


CHAPTER  Vllt  Dislocations  Moving  Uniformly  on  the  Interface 
Between  Two  Isotropic  Media  of  Different 
Elastic  PTOperties. 

The  problem  of  a  dislocation  moving  on  the  interface  separating 
two  media  of  different  elastic  properties  is  Interesting  both  from 
the  theoretical  as  well  as  the  practical  vleiypolnt.  Diffuslonlesa 
transformations  In  crystals  probably  Involve  dislocations  running 
on  the  Interfaces  between  transformed  and  untransfomed  material. 

Since  the  amount  of  energy  released  In  such  transformations  may  be 
large,  hl^  dislocation  velocities  are  to  be  expected.  In  fact, 
Eshelby  (1956)  has  proposed  that  dislocations  nay  run  at  supersonic 
velocities  In  diffusionless  transformations. 

A  dislocation  running  on  a  transformation  Interface  Is  moving 
on  an  Interfacd  which  separates  two  materials  of  differing  elastic 
properties  and  densities.  It  seems  worthwhile  to  bbtaln  the  solution 
of  the  stress  field  about  such  moving  dislocations.  This  section 
will  attempt  to  solve  the  problem  for  the  simplest  case:  that  in 
vdilch  the  two  elastic  media  are  isotropic.  (The  assusg>tlon  of 
isotropy  precludes  a  treatment  of  the  twinning  dislocation.  In 
isotropic  materials  a  twinning  dislocation  Is  merely  an  ordinary 
dislocation.  However,  the  analysis  may  have  some  qualitative 
application  to  twinning  dislocations  In  anisotropic  material.) 

Because  the  supersonic  velocity  range  may  be  of  practical  Importance 
the  dislocation  behavior  In  this  velocity  region  will  be  considered 
qvialltatlvely.  Consideration  Is  based  on  the  elucidating  analyses 


kh 


of  supersonic  dislocations  bjr  Esheltay  (1956)  cmd  Stroh  (1962) 
THEORY 


The  soluftion  of  the  problem  being  considered  can  be  obtained 
extending  the  known  solutions  of  dislocations  moving  In  Isotropic 
material  (Eshelby,  1949}  Frank,  1949;  Lelbfrled  and  Dietz,  1949  )• 


Consider  a  coordinate  system  in  which  a  dislocation  line  lies 
parallel  to  the  x  axis  and  moves  on  its  slip  plane  in  the  x  direction* 
Let  and  represent  the  Lame"'  constants  and  the  density 

of  the  material  above  the  slip  plane  (y  >  O),  and  ^2*  ^2 

the  same  constants  for  materials  below  the  slip  plane  (y  <  O).  In 
a  moving  screw  dislocation  the  following  equations  of  dynamic 
equilibrium  have  to  be  satisfied: 

1 


.2 
a  w . 


d  w 


1  _ 


„2 
a  w. 


(v.i) 


ax 


8t‘ 


oy  Ci 

where  1  =  1  or  2,  w^^  and  W2  are  the  elastic  displacements  in  the  z 
direction  above  and  below  the  slip  plane,  =  (p^/pj^)  '  -  the 

transverse  soimd  velocity  In  each  of  the  two  media,  and  t  stands 
for  time. 

A  moving  edge  dislocation  involves  the  following  equations 
of  dynamic  equilibrium: 


a  u 


a2 

a  u. 


^2 

a  V, 


a2 

a  u. 


V/  VZ  •  Vi  •  W  V  J  W  !*• 

+  2)ij^)  _____  +  _____  +  (kj^+  pj^)  _____  =  0^ 


(V.2) 


dx 


a2 

a  V 


ay 


a  V. 


5x8y 

a2 

a  u. 


+  2p^)  +  p.  +  -hxi) - =  Pl¬ 


at 

a 


ay 


ax 


Sxdy 


at* 
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idieve  u.  and  ▼.  are  the  elastic  displacements  In  the  x  and  7 
11  2 

directions  respectively.  The  tine  differential  __  can  be  replaced 

at^ 

2 

hy  for  the  dislocations  novlng  tri.th  a  uniform  velocity  V  in. 

0.2 

dx 

the  X  direction. 

From  the  results  of  analyses  of  dislocations  moving  in 
isotropic  material  ue  anticipate  that  the  solution  of  Bqns.  (V.l) 
and  (V.2)  for  moving  dislocations  will  bet 
(a)  Screw  Dislocation 


V  s 
1 


b 

2n 


tan 


-1 


(V.3) 


(b)  Bdge  Dislocation 
b 


"i  = 


^i  = 


2tt 

b 

4  n 


-1  Yiy  -1  ^ly 

tan  ^  tan  ^ 


^  log  (x'"^  +  Yi^y^)  +  log(x'^+ 


+  F. 


1 

*2. 


(V.4a) 


2  2^ 


In  these  equations  b  is  the  length  of  the  Burgers  vector, 
etc.,  are  constants,  \diere 

‘\i  =  t^^i'"^i^/Pl]  =  the  longitudinal  sound  velocity,  and 

=  x-Vt  where  V  is  the  velocity  of  the  dislocation.  (The  constant 
is  cuided  merely  to  make  Eqn.  (V.4b)  dlmensionslly  correct.  The 
constants  and  F^  are  added  in  order  to  match  suitably  the 
elastic  displacements  w^  and  u^  across  the  interface  y  s  0.  Since 
only  the  differentials  of  the  displacements  are  important  the 


(T-4b) 
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constants  can  be  Ignored). 

One  should  note  that  at  velocities  such  that  (or  y^)  i* 
an  inaglnary  number  the  arc  tan  of  (or  y^y/x  )  can  also  be 

viritten  as  log|(x'+  |frj^|y)/(x'  -  |Pj^|y)|  +  a  constant. 

The  evaluation  of  the  constants  for  the  screw  and  for  the  edge 
dislocations  are  considered  separately  tn  the  following  sections. 
SChlllU  DISLOCATIONS 

The  constants  of  equation  (V.3)  are  sinqple  to  evaluate. 

If  a  complete  circuit  is  made  around  a  dislocation  the  elastic 
displacement  must  change  ly  an  amount  equal  to  the  Burgers  vector. 
Thus  from  the  properties  of  the  arc  tan  function  it  is  evident  that 
+  Gg  =  2  when  the  dislocation  velocity  V  is  less  than  either 
c^  or  c^.  Another  equation  in  G^  and  G^  can  be  obtained  from  the 
condition  that  the  value  of  the  stress  must  be  continuous  across 


the  slip  plane.  Thus  at  z  =  0  the  stress  ™^Bt  satisfy  the 
condition  (<Jy2)i  “  ^'’y2^2* 


■  '“i 


A,  •‘1*’  vA 


/2  ^  .  2  2 

X  +  Pi  y 


(V.6) 


Therefore 

The  only  other  stress  which  exists  around  the  screw,  namely  d^,  is 
given  by  I 


-  .  i2  a  2  ( 

21T  X  +  Pi  y 


(V.8) 
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and  l8  equal  to  zero  on  the  slip  plane  regardless  of  the  values  of 
and 

A  solution  of  equations  (V.$)  and  (V.7)  results  int 


^2^2 

'l  "  - 

2PiPi 

-2  =  - 

*^1^1  >^2^2 


(V.9a) 


(V.9b) 


Equations  (V.9)  and  (V.3)  give  the  solution  of  the  elastic 
displacement  field  about  the  moving  screw  dislocation.  From  this 
field  both  the  stresses  (eqns.  (V.6)  and  (V.8))and  the  displacement 
velocities  can  be  found.  Once  these  quantities  are  known  the 
strain  energy  and  the  kinetic  energy  can  be  calculated  in  the  usual 
manner  (Kahelby,  1949;  Frank,  1949;  Leibfrled  and  Dietz,  1949; 
Weeirtman,  1961).  This  calculation  gives  for  the  total  self -energy 
E  of  a  screw  dislocation  moving  on  the  interface  between  two  laadla 


each  of  width  R  in  the  y  direction  and  extending  to  ±  »  in  the  x  direction 


E 


‘‘l‘‘2 


«  2  .  2^ 
*^2^2  ‘‘l^l 


3 


(V.io) 


As  would  be  e3q;>eoted  this  esqpresslon  becomes  infinite  at  the  slowest 
of  the  two  transverse  sound  velocities  (l.e.  when  or  P2  equal 
to  zero.) 
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SUPERSONIC  SCREW  DISLOCATION 


Two-dimensional  supersonic  dislocation  were  treated  first 
by  Eshelby  (1956)*.  Stroh  (1962)  later  more  generally  treated 
supersonic  dislocations  moving  in  anisotropic  media.  The 
difficulty  encountered  in  analyzing  supersonic  dislocations  is 
the  occvirrence,  when  linear  elasticity  theory  is  used,  of  Infinities 
In  the  equations.  This  problem  occurs,  for  example.  In  eqns.  (V.6) 
and  (V.8)  at  a  velocity  so  large  that  is  an  imaginary  number. 
According  to  these  equations  the  stresses  are  infinite  along  the 
planes  x*  =  y  •  Stroh  (as  well  as  Thomson  (1961))  points 
out  that  these  infinities  actually  do  not  occur.  Because  of  the 
discrete  atomic  nature  of  a  crystal  lattice,  only  finite  stresses 
and  finite  energies  can  exist.  It  could  be  assvoaed  that  linear 
elasticity  theory  holds  only  for  those  stresses  ’■diose  "effective 
shear  stress"  t  (  defined,  for  example,  by  Nye  (1957)  as  the  sqvutre 

-  is  the  Kroneker 


root  of 


2  ^ 
^  ij 


where  a 


I  _ 


ij  '^Ij 


delta,  and  is  the  usvial  stress  component)  does  not  exceed  the 

value  X  .  Screw  dislocations  when  only  the  stresses  c  and  a 

o  ^  yz  xz 

are  present  may  be  expressed 

C  ^Z  +  d  ^z  +  '"o^ 
y  X 

The  breakdown  in  linear  elasticity  theory  in  supersonic  situations 
occurs  along  the  planes  x'  =  ±  |  p^  |  y  when  p^  becomes  imaginary. 
These  discontinuities  in  the  displacements  and  stresses  are  the 
plane  waves  vdiich  Eshelby  and  Stroh  have  shown  must  exist  in  order 
to  maintain  a  supersonic  dislocation. 

Figure  V.l  shows  schematically  the  behavior  of  the  Eshelby-Stroh 


*0ne-dimen8ional  supersonic  dislocations  were  studied  earlier  by 
Frank  and  van  der  Merwe  (1950) 
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discontinuity  planes  as  the  velocity  of  a  screw  dislocation  is 
increased.  The  lines  in  this  figure  represent  the  region  \diere  the 
stress  equals  the  maxlimun  stress  For  velocities  less  than 
the  stress  is  equal  to  only  in  the  core  region.  At  the  velocity 
0^  the  Eshelhy-Stroh  discontinuity  appears.  The  normal  to  this 
discontinuity,  which  moves  with  the  velocity  c^,  is  in  the  direction 
of  the  dislocation  and  hence  the  discontinuity  dissipates  no 
energy  to  the  crystal  svirface.  At  a  somevdiat  greater  velocity 
two  K-S  discontinuities  appear.  The  directions  of  their  normals 
are  different  from  the  direction  of  dislocation  motion  and  hence 
surface  tractions  are  reqvdred  to  maintain  the  discontinuities. 

If  svdtable  surface  tractions  are  not  applied,  the  dislocation 
still  can  run  but  the  slip  plane  must  be  able  to  give  up  energjr. 

It  could  do  so  either  if  an  external  shear  stress  is  applied  to 
the  crystal  or  if  the  slip  plane  is  slso  a  transformation  plane 
which  gives  up  energy  as  the  dislocation  runs  along  it.  Vftien  no 
surface  tractions  are  applied  the  E-S  discontinuities  which  are 
ahead  of  the  discontinuity  behind  the  dislocations  will  occur. 

This  E-S  discontinuity  sends  energy  out  to  the  crystal  surface 
and  an  equal  amount  of  energy  must  be  supplied  at  the  dislocation 
core  if  the  dislocation  is  to  continue  to  run. 

It  is  simple  to  obtain  an  estimate  of  this  energy 
dissipation.  The  E-S  discontinuity,  which  is  expected  to  leave 
a  width  of  the  order  of  the  atomic  spacing,  has  associated  with 
it  an  energy  per  imlt  area  of  the  order  of  x^b.  The  discontinuity  makes  an 
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langle  with  the  surface  of  the  order  of  tan~^  |l/Pj^|.  When  a 
unit  length  of  dislocation  moves  a  unit  distance, an  amount  of 
energy  equal  to  Is  dissipated  and  thus  the  E->S  discontinuity 

produces  a  retarding  stress  c  equal  to 

(V.12) 

Since  r^  Is  the  theoretical  shear  stress.  It  seems  to  be  Impossible 
to  get  to  stresses  much  above  c^.  At  greater  velocities  the 
retarding  stress  Is  greater  than  the  stress  which  the  slip  plane 
can  support* 

Because  of  the  retai^ilng  stress  o^,  it  may  not  be  possible 
to  bring  a  dislocation  from  zero  velocity  gradually  up  to  and 
beyond  the  velocity  c^.  If  the  velocity  c^  Is  much  smaller  than 
c^*  Equation  (V.6)  shows,  however,  that  If  the  dislocatlion  Is 
set  running  at  the  velocity  c^,  the  displacements  above  the  slip 
plane  can  be  set  equal  to  zero.  An  discontinuity  thus  exists 
only  below  the  slip  plane,  as  shown  In  Flg.V.l.  Since  the  normal 
of  this  discontinuity  is  in  the  direction  of  dislocation  motion, 
no  energy  dissipation  is  involved  and  the  dislocation  can  run 
without  energy  dissipation.  (Other  energy  dissipation  mechanisms, 
such  as  that  due  to  dispersion  (Eshelby,  19$6,)  will  still  operate.) 
The  velocity  c^  thus  is  singular  in  that  the  dislocation  again  acts 
like  a  subsonic  dislocation.  This  type  of  singular  dislocation 
velocity  in  the  supersonic  region  was  first  noticed  by  Eshelby  (1949) 
for  edge  dislocation.  This  type  will  be  considered  again  in 
the  next  section  on  edge  dislocations.  Eshelby 's  singular 
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dislocation  velocity  probably  approximates  the  velocity  at  \Aiioh 
diffusionless  transformations  take  place.  It  represents  a  fast 
dislocation  velocity  at  vdiich  energy  dissipative  processes  from 
the  E-S  discontinuities  are  minimized.  (Of  course  other  dissipative 
processes  still  can  occur.) 

Above  the  velocity  C2  the  E-S  discontinuities  reappecur  above 
the  slip  plane,  as  sho\m  in  Fig.  V.l. 

EDGE  DISLOCATION 

The  problem  of  an  edge  dislocation  moving  on  an  interface 
separating  tw  different  elastic  media  is  more  conqpllcated  but 
more  interesting  than  that  of  the  moving  screw  dislocation.  The 
constants  appearing  in  Eqns.  (V.4)  for  the  elastic  displacements 
are  slnqply  evaluated.  If  Eqns.  (V.4)  are  placed  into  Eqns.  (V.2) 
one  finds  that  (V.13a) 

\  =  Pi"\  (V.13b) 

The  condition  that  the  displacements  describe  an  edge  dislocation 
gives  the  equation 

JIj^  +  Bj^  +  A2  +  B2  =  2  (V.14a) 

The  condition  that  no  line  forces  act  at  the  dislocation  core  results 
in 

‘"l^Vl^  +  Bj^)  +  1^2^  V2^  V  "  ° 

where 

The  requirement  that  the  stresses  be  continuous  across  the  slip 
plane  produces  the  equation 
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Finally I  the  condition  that  the  elaetic  displacements  be  contizxuous 
across  the  slip  plane  gives 


Hhen  Eqns.  (V.I4)  are  solved  for  etc.,  oxie  obtains 


AP,P 


1*^2 


“l^»^2^l^l“l  "  ^2^2  ^  ^2*^2  ^>^1^1 


•*•  >^2^2  “  t^lPl^l^  “  ^lh'f2^'^l 


(V.15a) 


h  = - }  '(\^2^2^^2^2  -  >^1^  +  W2P2<»^I  "  1^2^ 


AP, 


+  -  1)) 


(V.15b) 


vbere 


A  = 


PiP 


-  (p2''2^  “  »*l“l  -  ^‘l^*2^Yl  ^  Y2^^Pi‘‘2  ■*■  P2“l  ^ 


1*'2 


+  (p|^ri  +  ■*■  *^2^2^  ”  YiYaPiPg^i^i  " 

-  ^^l»‘2<Pl  +  p2^<Yi''2^  +  Y2^i^)  + 

To  obtain  A^  and  one  merely  interchanges  the  subscripts 
1  and  2  in  these  equations. 

STRESS  ON  THE  DISLOCATION  SLIP  PLANE  AND  THRESHOLD  DISLOCATIONS  VELOCITY 
With  these  values  of  the  constants  A^,  etc.,  one  obtiixis  the 
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following  expression  for  the  shear  stress  o  aotlng  on  the  slip 
plane. 


(y,(1 

pc  PiPai  ( 


2 

“2 


)jii(yiPi 


+ 


Yi(l-ai  )P2^Y2P2  " 


(V.16) 


It  has  been  found  from  a  study  of  edge  dislocations  moving  In 
an  Isotropic  material  (Ueertmanj  1961)  that  as  the  dislocation 
velocity  Is  increased  the  shear  stress  on  the  slip  plane  decreases 
until  It  becomes  zero  at  the  Raylel^  wave  velocity.  At  greater 
velocities  It  increases  with  Increasing  velocity  but  has  a 
negative  value.  For  the  type  of  dislocation  under  considerations 
the  velocity  at  vdiich  the  shear  stress  on  the  slip  plane  goes 
to  zero  can  be  found  by  setting  Fqn.  (V.16)  equal  to  zero.  The 
following  equation  Is  obtained: 


(YiPi  -  _  (Y3P2  -  “2^^  =  ° 

YiPi  Y2P2 

The  velocity  which  satisfies  that  eqviatlon  Is  the  threshold 
velocity  (Teutonlco's  term  (1962b))  separating  the  region  of 
normal  dislocation  behavior  (in  which  dislocations  of  like  sign 
on  the  sane  slip  plane  repel  one  another)  from  the  region  of 
anomalous  behavior  (dislocations  of  like  sign  attract  one  another). 

The  equation  Yj^Pj^  “  ~  0  (V.18) 

determines  the  Rayleigh  surface  wave  velocity  In  the  Isotropic 
medium  above  or  below  the  slip  plane.  Henoe  one  can  see  from 
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Eqn>  (V.17)  that  the  threshold  velocity  lies  between  the  Rayleic^ 
wave  velocities  of  each  medliim*  The  threshold  velocity  defined 
by  Eqn.  (V,17)  is  not  equal  to  the  Stoneley  wave  velocity*,  contrary 
to  idiat  one  might  have  assmed  from  the  fact  that  the  threshold 
velocity  is  equal  to  the  Rayleigh  wave  velocity  when  the  elastic 
properties  and  densities-  of  the  two  media  are  identical. 

It  is  not  always  possible  to  find  a  velocity  which  satisfies 
Eqn.  (V.17) .  Figure  V.2  Illustrates  this  situation  for  the  simple 
case  in  \dilch  the  longitudinal  sound  velocity  in  each  medlmn  is 
very  much  larger  than  both  Cj^  and  ^2  so  that  both  Yg  are 

equal  to  1  for  velocities  near  c^  axid  Cg.  The  densities  in  the 
two  media  are  assiuaed  to  be  eqvial.  One  can  see  from  the  figure 
that  the  ocoxirrence  of  a  threshold  velocity  is  limited  to  a  narrow 
range  in  the  variables  and  pg  varying  from  1  to 

1»193)»  If  ^8  outside  this  range  the  anomalous  velocity 

region  cannot  exist. 

SUPERSONIC  EDGE  DISLOCATION 

Above  the  slower  of  the  two  transverse  sound  velocities  the 

edge  dislocation  is  in  a  supersonic  region.  As  in  the  ease  of  a 

screw  dislocation  Eshelby-Stroh  discontinuities  will  arise  along 

the  planes  x'  =  ±|  1 7  vrtien  Yj^  also  is  an  imaginary  quantity. 

The  behavior  to  be  ejected  is  shown  in  Fig.  V.3  for  the  case  in 

which  There  are  no  discontinuities  below  the 

*Stonely  waves  (Ewing,  Jardetzky,  and  Press,  1957)  are  the  surface 
waves  vAilch  propagate  along  an  interface  separating  two  eeml-lnflnlte 
media  of  differing  elastic  properltes.  They  are  slsqply  a  generalisation 
of  the  Rayleigh  wave  vdiich  propagates  along  a  free  surface. 
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X/2 

FIG,  V.2  Plot  of  the  threahold  velocity  versus 

The  velocity  Is  expressed  In  units  of  o^*  Also 
shoun  Is  the  variation  of  the  transverse  sound 
velocities  c^  and  c^  and  the  two  Rayleigh  surface 
wave  velocities  c^^^  and  c^^. 
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velocity  except  at  the  dislocation  core.  At  the  velocity 
and  E>^  discontinuity  occurs  vdiose  normal  is  the  direction  of 
dislocation  notion.  In  contrast  to  the  case  of  the  screw 
dislocation,  the  stresses  about  an  edge  dislocation  in  the  upper 
plane  are  (when  linear  elasticity  is  applied)  infinite  eveiy- 
idiere  rather  than  in  a  region  limited  to  the  plaine  x'  =0. 

This  infinity  in  energy  is  thus  of  a  different  kind  from  that 
of  a  screw  dislocation.  Stroh  has  pointed  out  that  the  difference 
is  due  to  a  resonance  phenomenon.  At  velocities  between  c^  and  c  ^ 
the  normals  of  B-S  discontinuities  make  an  angle  with  the  direction 
of  dislocation  motion.  The  dislocation  experiences  a  retarding 
•force  given  again  by  Bqn.  (V.12).  At  higher  velocities  other  E-S 
discontinuities  appear,  as  shown  in  the  figure. 

ESHELBY'S  SINGULAR  DISLOCATION  VELOCITY 


Eshelby  (1949}  showed  that  an  edge  dislocation  in  an 
isotropic  medium  can  move  at  a  velocity  J  2c  without  any  E-8 
discontinuity  appearing  in  its  displacement  field.  This  velocity 
is  in  a  supersonic  range  and  yet  the  dislocation  will  not  radiate 
energy.  The  reason  why  the  dislocation  can  exhibit  behavior  can 
be  seen  from  Bqns.  (V.U;.)  for  the  case  of  ”  ^2* 

=  pg.  In  this  situation  Eqns.  (V.14c)  and  (V.14d)  are  identically 
aero  since  A^  =  and  Now  and  of  Eqns.  (V.14a)  and 

(V.14b)  could  be  set  equal  to  zero  (and  thus  eliminate  the  terms 
of  Eqn.  (V.4)  which  contain  the  E»S  discontinuity  vdien  0  <  V  <  0^) 
and  still  these  equations  may  be  satisfied  provided  the  velocity 
such  that  a^^  “  ^2  ”  velocity  is  ,^20.  At  this  velocity 
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FIG.  V.3.  Sohematlo  plot  of  the  Bahelbgr-^troh  dlsoontlxmitles 
for  a  moving  edge  dislocation  at  different 
dislocation  velocities.  It  is  assumed  that 

°1  ®2  ‘Ovl  ®K2* 


Qo  discontixxultieB  appear  in  the  stress  and  displacement  field 
and  all  of  Eqn.  (V.14)  are  satisfied. 

Consider  now  vdiat  happens  to  Eshdlbgr's  singular  dislocation 
velocity  in  the  edge  dislocation  being  discussed.  For  this 
singular  velocity  to  esdst,  it  must  be  possible  to  set  any  of  the 
and  of  Eqns.  (V.14)  equal  to  zero  when  the  terns  multiplied 
by  axid  of  Eqns.  (V.4)  contain  E-S  discontinuities.  If  the 
dislocation  velocity  lies  between  c^  and  c^i  and  B^,  therefore, 
is  set  equal  to  zero,  Eqns.  (V.  4)  caimot,  in  general,  hold 
for  any  vtdue  of  V  since  we  have  four  equations  in  three  unknown 
qiiantitles.  However,  it  is  possible  for  all  equations  to  hold 
if  the  determinate  of  the  coefficients  of  Eqns.  (V.14b)  throu^ 
(V«14d)  is  zero.  Thus  the  following  equation  for  Eshelby's 
singular  dislocation  velocity  is: 

—  ^'^2^2  “  “2^^  “  ^^2^2  “  “2^^  "  “2^^  “  ° 
lil  Yi 

For  the  case  in  \dilch  ”  Pg  very  much  larger 

than  both  c^^  and  Cg  so  that  as  =  Y2  ~  ^  velocities 

of  interest,  it  is  found  that  as  varied  from  1  to  infinity 

the  Eshell^  singular  velocity  c_  varies  as  shown  in  Fig.  V.4* 

<9 

The  singular  velocity  can  exist  when  '  P2*  ^2  sUfi^tly 

larger  than  no  singular  velocity  occurs.  In  the  range 
2  ^  ^  4  *  singular  velocity  does  exist.  Its  value  is  almost 

equal  to  the  velocity  o^.  In  contrast  with  the  case  df  the 
Eshelby  singularity  in  the  screw  dislocation,  this  singular  velocity 
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units  of 


has  no  Infinite  stresses  associated  with  it. 

If  =  1^2  i'*'  1-8  not  possible  to  have  a  singular  velocity 
vMch  is  larger  than  both  c^  and  C2*  However,  if  =  |ii2 

is  not  necessarily  equal  to  P2)  it  is  possible  to  have  a  singular 
velocity  greater  than  c^  and  C2.  In  Eqns.  (V.I4)  and  B2  can 
be  set  equal  to  zero.  V/hen  this  is  done,  one  obtains 


(V.20) 

For  Eqn.  (V.14b)  to  hold  the  velocity  must  be  such  that 
a/+  (Yi/Y2)a/  =  0  (V.2l) 

This  equation  defines  the  singular  velocity.  When  =  Y2  ~  ^ 
this  velocity  is  given 


2  _ 


,22  , 

4Cj^  C2  4^1 


(V.22) 


«1  ^«2 


Pl  +  P, 


SUI4MARX 

The  behavior  of  dislocations  moving  on  the  Interface  between 
two  different  elastic  media  generally  is  fomid  to  be  what  one 
would  expect  from  the  studies  of  dislocations  moving  in  isotropic 
material.  The  slowest  sound  velocity  sets  an  upper  limit  (apart 
from  the  Eshelby  singular  velocity)  to  the  speed  of  the  dislocation 
since,  according  to  linear  elasticity  theory,  the  energy  becomes 
infinite  at  this  point.  If  it  is  assuad  that  a  crystalline  material 
can  support  only  a  finite  stress,  it  can  be  shown  that  supersonic 
solutions  exist.  Eshelby  and  Stroh  have  pointed  out  this  fact. 
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Because  of  damping  caused  by  the  geneiratlon  of  sound  uaves»  the 
supersonic  dislocation  velocities  probably  are  limited  to  Teloclties 
not  much  beyond  times  a  sound  velocity. 

The  Eshelby  singular  dislocation  velocity  may  or  may  not 
exist.  Its  existence  depends  on  the  values  of  the  elastic  ooxietants 
and  the  density  in  each  medium.  Because  energy  dissipation  processes 
are  minimized  at  the  Eshelby  singular  velocity,  this  velocity  could 
be  that  at  which  transformation  dislocations  run  on  diffusionless 
transformation  interfaces.  It  would  be  Interesting,  therefore,  to 
see  if  fast  diffusionless  transformations  take  place  in  materials 
whose  elastic  constants  do  not  permit  the  existence  of  an  Eshelby 
singular  velocity. 

The  threshold  velocity  separating  a  nonnal  from  an  anomalous 
velocity  region  also  may  or  may  not  exist.  It  will  exist  only 
if  the  two  transverse  sound  velocities  In  the  two  isotropic  media 
above  and  below  the  slip  plane  lie  close  to  each  other. 
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AFB,  Ohio 

(ASAPRL,  Technical  Doc  Library) 

(ASRMDS-l,  Mr.  Janik) 
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BSD  (Ballistic  Systems  Division),  Norton  AFB,  Calif 
(Tech  Library) 

(BSR) 

(BSRD,  Lt  Col  Caseria) 

(BSRVE,  Lt  Col  Parker) 

(BSLA) 

(BSTD) 

(BSAT) 

SSD  (Space  Systems  Division),  (SSSC-TDC),  AF  Unit  Post  Office, 
Los  Angeles  45,  Calif 

ESD,  Hanscom  Fid,  Bedford,  Mass 

(ESAT) 

(ESDL,  Col  R.  J.  Lynch) 

(ESDS,  Col  W.  H,  Congdon) 

AF  Msl  Dev  Cen,  (AFOSR,  SRAT),  Holloman  AFB,  NM 

AFFTC  (FTFT),  Edwards  AFB,  Calif 

AFMTC,  (Tech  Library,  MU-135),  Patrick  AFB,  Fla 

APGC  (PGAPI),  EglinAFB,  Fla 

RADC  (Document  Library),  Griffiss  AFB,  NY 

KIRTLAND  AFB  ORGANIZATIONS 
AFSWC,  KirtlandAFB,  NM 
(SWEH) 

(SWOI) 

(SWRPA) 

(SWRPL) 

(SWRPT) 

(SWRA) 

ADC  (ADSWO),  Special  Weapons  Office,  KirtlandAFB,  NM 

ATC  Res  Rep(SWN),  AFSWC,  KirtlandAFB,  NM 

AFLC,  Albuquerque  Ln  Ofc  (MCSWQ),  AFSWC,  KirtlandAFB,  NM 

SAC  Res  Rep  (SWL),  AFSWC,  KirtlandAFB,  NM 

TAC  Liaison  Office  (TACLO-S),  AFSWC,  KirtlandAFB,  NM 

US  Naval  Weapons  Evaluation  Facility  (NWEF)  (Code  404), 

Kirtland  AFB,  NM 
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OTHER  AIR  FORCE  AGENCIES 

Director,  USAF  Project  RAND,  via:  Air  Force  Liaison  Office, 

The  RAND  Corporation,  1700  Main  Street,  Santa  Monica,  Calif 

{RAND  Library) 

Physics  Dept  (Dr,  Olen  Nance) 

Physics  Dept  (Capt  J,  A,  Welch) 

Physics  Dept  (Mr,  Jack  Whitener) 

Aerospace  Defense  Systems  Office  (ADO),  ATTN;  ADSO, 

Air  Force  Unit  Post  Office,  Los  Angeles,  Calif 

ARMY  ACTIVITIES 

Chief  of  Research  and  Development,  Department  of  the  Army 
(Special  Weapons  and  Air  Defense  Division),  ATTN:  MaJ  Baker, 

Wash  25,  DC 

Commanding  Officer,  Diamond  Ordnance  Fuze  Laboratories, 

(ORDTL  06.33,  Technical  Library),  Wash  25,  DC 

Director,  Office  of  Special  Weapons  Development  (USACDC), 

Ft  Bliss,  Tex 

ARGMA  Liaison  Office,  Bell  Telephone  Labs,  Whippany,  NJ 

Redstone  Scientific  Information  Center,  US  Army  Ordnance  Missile 
Command  (Tech  Library),  Redstone  Arsenal,  Ala 

Director,  Ballistic  Research  Laboratories,  Aberdeen, 

Proving  Ground,  Md 

Library  (Mr,  Ed  Bailey) 

Tenninal  Ballistics  Lab  (Dr,  Coy  Glass) 

Commanding  Officer,  Picatinny  Arsenal,  Samuel  Feltman 
Ammunition  Laboratories,  Dover,  NJ 

(SMVPA-VA6) 

OHD  BBVC2  (Mr,  Murray  Weinstein) 

Research  Analysis  Corp,  (Document  Control  Office),  6935  Arlington 
Road,  Bethesda,  Md,  Wash  14,  DC 

US  Army  Office  of  Ordnance  Research  (Physical  Sciences  Division), 
Box  CM,  Duke  Station,  Durham,  NC 

Director,  Army  Research  Office,  Arlington  Hall  Sta,  Arlington,  Va 

Commanding  Officer,  US  Army  Engineers,  Research  &  Development 
Laboratories,  Ft  Belvoir,  Va 

Commanding  General,  White  Sands  Missile  Range  (Technical  Library), 
White  Sands,  NM 
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NAVY  ACTIVITIES 

I  Chief  of  Naval  Operations,  Department  of  the  Navy,  (OP-36), 

Wash  25,  DC 

I  Chief  of  Naval  Research,  Department  of  the  Navy,  (Mr.  James 

Winchester,  Wash  25,  DC 

1  Chief,  Bureau  of  Naval  Weapons,  Department  of  the  Navy, 

Wash  25,  DC 

1  Commanding  Officer,  Naval  Research  Laboratory,  Wash  25,  DC 

2  Director,  Special  Projects,  Department  of  the  Navy,  (Mr.  D.  R. 
Williams,  Wash  25,  DC 

OTHER  DOD  ACTIVITIES 

Chief,  Defense  Atomic  Support  Agency,  Wash  25,  DC 
2  (Document  Library) 

1  (DASARA,  Lt  Col  Singer) 

1  Commander,  Field  Command,  Defense  Atomic  Support  Agency 

(FCAG3,  Special  Weapons  Publication  Distribution),  Sandia  Base, 
NM 

1  Director,  Weapon  Systems  Evaluation  Group,  Room  2E 1006, 

The  Pentagon,  Wash  25,  DC 

5  A^va^c^dj Research  Projects  Agency,  The  Pentagon,  Wash  25, 

1  Director,  Defense  Research  &  Engineering,  (Col  Gilbert),  The 

Pentagon,  Wash  25,  DC 

10  ASTIA  (TIPDR),  Arlington  Hall  Sta,  Arlington  12,  Va 

1  DASA  Data  Center,  Tempo-General  E  lectric  Co.,  P.  O.  Drawer 

QQ,  Santa  Barbara,  Calif 

AEC  ACTIVITIES 

1  US  Atomic  Energy  Commission,  (Headquarters  Library), 

Wash  25,  DC 

2  Sandia  Corporation,  (Technical  Library,  Mrs.  Allen),  Sandia 
Base,  NM 

2  Sandia  Corporation,  (Tech  Library,  C.  A.  Paxton),  P.  O.  Box  969, 

Livermore,  Calif 

1  Chief,  Division  of  Technical  Information  Extension,  US  Atomic 

Energy  Commission,  Box  62,  Oak  Ridge,  Tenn 

1  University  of  California  Lawrence  Radiation  Laboratory  (Technical 

Information  Division,  ATTN:  Mr.  Clovis  Craig),  P.  O.  Box  808, 
Livermore,  Calif 
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University  of  California  Lawrence  Radiation  Laboratory,  Technical 
Information  Division,  Berkeley  4,  Calif 

Director,  Los  Alamos  Scientific  Laboratory,  (Helen  Redman. 
Reports  Library),  P.  O.  Box  1663,  Los  Alamos,  NM 

Brookhaven  National  Laboratory,  Upton,  Long  Island,  NY 

Argonne  National  Laboratory  (Tech  Library),  Argonne,  Ill 

Oak  Ridge  National  Laboratory  (Tech  Library),  Oak  Ridge,  Tenn 

OTHER 

Battelle  Memorial  Institute,  505  King  Avenue,  Columbus,  Ohio 

Institute  for  Defense  Analysis,  Room  2B257,  The  Pentagon, 

Wash  25.  DC 

Institute  of  the  Aerospace  Sciences,  Inc. ,  2  East  64th  Street, 

New  York  21,  NY 

Colorado  School  of  Mines,  Golden,  Colo 

ATTN:  Dr.  John  Rinehart  and  Capt  R.  A.  Houidobre 

Republic  Aviation  Laboratory,  MSVD,  (Dr.  Tcutonico),  223  Jericho 
Turnpike,  Mineola,  Long  Island,  NY 

Aerojet-General,  (Mr.  Waldorf),  Downey,  Calif 

General  Electric,  Space  Sciences  Laboratory,  (Dr.  D.  Flom), 
Valley  Forge.  Pa 

Stanford  Research  Institute,  Poulter  Laboratory,  (Dr.  George 
Abrahamson),  Menlo  Park,  Calif 


* 


> 


1  The  Boeing  Company,  Aero-Space  Division,  (Mr,  Penning), 

P.  O.  Box  3707,  Seattle  24,  Wash 

1  Armour  Research  Foundation,  (Mr.  Dennen),  10  West  25th  Street, 

Chicago,  111 

1  Avco  Corporation,  201  Lowell  Street,  Wilmington,  Mass 

ATTN:  Dr.  D.  Morgan 

1  Official  Record  Copy,  (SWRPA,  Capt  Lamberson) 

1  OTS,  Department  of  Commerce,  Wash  25,  DC 
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